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This paper develops a hierarchically-structured, suboptitnal controllt^r 
for a linear stochastic system composed of fast snd slow subsystems. It is 
proved that the controller is optimal in the limit as the separation of 
time scales of the subsystems becomes infinite. The methodology is illus- 
trated by design of a controller to suppress the phugoid and short period 


modes of the longitudinal dynamics of the F-8 aircraft. 
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A common occurrence in engineering systems is tlie presence of 
phenomena that naturally evolve in widely separated time scales. So"i 
examples drawn from the fields of aerospace and power engineering w' 
serve to illustrate this point. 

It is well-known [1] that the longitudinal dynamics of an aircra.t 
are comprised of two distinct oscillatory modes - the phugoid and short 
period - of periods on the order of 100 sec. and 1 sec. A terrestrial in- 
ertial navigator has Schuler oscillations of period 84 min. and earth rate 
oscillations with period 24 hr. [2]. A dual spin satellite in synchronous 
orbit will be subject to an orbital oscillation with a 24 hr. period and 
a nutation with period on the order of 1-10 sec. [3,4], These effects are 
often used in an ad hoc way in the design of filters and controllers, 
usually by assuming that the slow modes are constant if the fast modes are 
of concern, or by ignoring the fast modes when the slow modes are of inter- 
est. See [3, 5-8, 28] for examples. 

Additional examples can be found in the field of electric power 
systems. An electrical machine has an oscillatory node involving stator 
fluxes that is invariably neglected in favor of the much slower electro- 
mechanical oscillations in stability studies [9] . A similar approximation 
is made in studies of a large number of interconnected machines, in which 
the intermachine electrcmechanical swings are ignored when the much slower 
average frequency behavior is of primary concern [10] . 

In addition to these concrete examples, note that proponents of 
hierarchical control often suggest that the task of controlling a large 
scale system should be partitioned into sxdetasks by tim.e scale. Thus the 
higher levels of the control system are concerned with slower phenomena, 
and the lower levels with faster phenomena [11-14]. It is difficult to 
point to any specific examples, with the possible exception of the inter- 
action’ between automatic generation control and economic dispatch on 
Joctric power systems [15]. 

The multiple time scale phenomena alluded to above are conveniently 
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modelled via perturbation theory (16], There are a number of possible ap- 
proaches, but we will adopt the framework of singular perturbation theory. 
This theory has been applied to a variety of control problems by a number 
of authors (17-21, 24-27) but the previous work most relevant to this paper 
is that of Kokotovic et. al. [18-20] and Haddad (21). This paper is in the 
spirit of [18-21], and moreover requires several of the detailed results of 
these papers concerning singular perturbations of Ricatti equations. 

Specifically, the paper begins by analysis of singular perturbations 
for linear stochastic systems with two time scales. An approximating system 
is obtained with the property that the mean-square error between the states 
of the actual and approximating systems approaches zero as the separation 
of time scales becomes infinite. The usefulness of this result is demon- 
strated by application to the stochastic optimal linear regulator problem. 
With the machinery properly set up, it is straightforward to identify the 
asymptotically optimal controller, using known results on singular pertur- 
bations of Riccati equations. The controller has an Interesting hierarchical 
structure, with the implication of reduced on line computations. 

These new theoretical results are illustrated by application to an 
important control problem. An asymptotically optimal two time scale control- 
ler is developed for the longitudinal dyneunics of a jet aircraft. The two 
time scale controller is compared to the optimal controller, and it is 
demonstrated that there is negligible degradation in performance. 

An attempt is made throughout to relegate technical details to the Ap- 
pendix, so that the paper will be accessible to engineers interested only 
in the two time scale design procedure. 

MAIN RESULTS 

Singulau: perturbation theory is concerned with systems of the form 


x(t) e) ■ f(x, y, e) (2.1) 

e^(t; e) - g(x, y, e) (2.2) 

and the corresponding degenerate system 

x(ti 0) - f(x, y, 0) (2.3) 

0 ■ g(x, y, 0) (2.4) 

The basic question is whether (2.3) - (2.4) is an approximation to 
(2.1) - (2.2) In the sense that 

lim x(t; e) - x(t; 0) (2.5) 

£-►0 


lim y(t; e) ■ y(t; 0) (2.6) 

e -^0 

Various technical assumptions are required to obtain (2.5) emd (2.6), but 
under these assumptions the degenerate system is a valid reduced order ap- 
proximation to the original system in the sense that for E sufficiently 
small the solutions of the two systems are close. 

In the stochastic case, the situation is more complex. Consider the 
linear system 


dt 


x^(t; e) 
ex2(t» E) 


Aii(f:) X^(t; E)1 Lj^(e) 

I ♦ C(t) (2.7) 

A2i(E) A22(e)J [x2(t; E)J [^2^^^^- 


2 . 


where 


E 


Xj^(0> e) 

■ 

*2^0; C) 

T 

x^(0; e)Xj^(0; 
XjCO; e)x^(Oj 


0 

0 

m M 

T 

E) Xj^(Oj e)x2(0i E) 

e) x^'O; E)x2(0; E) 


z z 

^11 ^12 


z z 

L 21 22 


E{C(t)} - 0 
E{^(t)C'^(s)} - H6(t-s) 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 
( 2 . 11 ) 


Xj^(O), X 2 ( 0 ) aure independent of ^^(t) , matrices are continu- 

ous in E at £ * 0. Moreover, ^22^^^ stable. 

An approximation to (2.7) is desired that is valid for small C and is 
simpler than (2.7). Note that setting E ■ 0 in (2.7) is inadequate; since 

x^(tj 0) - - A22^(0)A2^(0)x^(t; 0) - A22^ (0) L 2 (0) C (t) (2.12) 


has a white noise component and therefore has infinite variance. Consequent- 

ly, 


X2(t; 


e) - 


X 2 (t; 


»)’( 


X 2 (t; 


e) - 


X 2 (t 


; 0 ) 


+ 00 


(2.13) 


so that X 2 (t; 0) is not am approximation to X 2 (t; e) (in the least squares 
sense) . 


Instead, define the stochastic degenerate system associated with (2.7) 
to be the system 


where 


A 

L 

A 

L 



e) 

■ *lld 

(E)Xj^^(t; 

E) + Lj^^(E)C(t), x^^(0; 

E) - x^(0) 

(2. 

.14) 

•0 

CM 

E) 

■ ^Id 

(E)Xj^^(t; 

E) + ^22d^^^*2d^^’ ^ 

L2^(E)C(t), 





■‘2d<°‘ • *2 

(0) 


(2. 

,15) 


- 

A^l(e) 


22^‘«*2l'« 


(2. 

.16) 

ld‘^> 

- 

Lj^(e) - 

Ai2(E)A22 

^(E)L2(E) 


(2. 

.17) 

21d*'^’ 

- 

A2i(C) 




(2. 

.18) 

2d'^> 

- 

Lj(e) 




(2, 

.19) 


Notice that the stochastic degenerate system is of the same order as the 
original system, unlDce the situation for deterministic singular pertur- 
bations. 


Theorem 1 


Consider the linear stochastic system (2.7) - (2.11) and a correspond- 


I 


ing stochastic degenerate system (2.14) - (2.15). Assume that all matrices 
in the two systems are continuous in e at e •* 0, and that A (0) and 
-1 1 

^ 11 ^®^ ” ^12 ^^^^22 ( 0 )^ 21 ^®^ stable (i.e., have eigenvalues in the 


open left half complex plane) . Then the stochastic degenerate system is an 
approximation to the original system in the sense that 


lim 

lim 

0 



e) 

- -Id''’ 

G) ) (Xj^ (t; 

e) 

- 

e))"^} - 0 

(2.22) 

E{(X2(t; 

e) 

- ’“2d"' 

E) ) (X2(t; 

e) 


E))"^} - 0 

(2.23) 


uniformly for 0 < t < ®. 


Proof 

The proof of the theorem is quite involved. Differential equations for 


Ell(t;e) - E{Xj^(t; e)x‘ (t; e) } - E{(Xj^(t; e)-x^^(t;E)) (*j(t|£)-Xj^^(t; e) )^} 

(2.24) 

222(t;£) ■ eCx^ (t; e)x2(t; e) } ■ E{ (x^ (t; e) -x^^ (t; £) ) (x^ (t; e) -x^^ (t; G) ) '^} 

(2.25) 

are obtained, and the limits are evaluated by (non-stochastic) singular 
perturbation theory to establish (2.22) and (2.23). See the Appendix for 
details. 

Remeurks 

1. Note that the stochastic degenerate system is the same order as the 
original system, so that (2.23) is valid for t 0. 

2. Clearly, the assunption ^^^^^(0) - stable is only 

necessary to insure uniform convergence in (2.22), (2.23) on the infinite 
interval. Without this assumption, a theorem of Tihonov (16) can be invoked 
which insures uniform convergence in (2.22), (2.23) for sets of the form 
I0,T). 

3. Proof is easily generalized to cover uniformly asymptotically 
steUsle time- vary ing systems at the expense of some additional notation. 

Consider now the system ^ 

Jc(t; e) - A(e)x(t; e) B(e)u(t) + L(e)C(t) (2.26) 

with observations 

y(t; e) “ Cx(t; E) + 0(t) (2.27) 

and cost 

r 

J - liin ~ IE{ I x^Qx + u*^Ru dt}] (2.26) 

T*- 00 O’' 


Note that this assumption implies that there exists an £q > 0 such that the 
systern matrix in (2.7) is stable for all 0 ^ E < £q (22). 
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where 


^12 

^21 ^22 


» B - 


1 ®21 

Cl C^J , Q. 

_ ®21 ®22 


EU(t)} » 0, 


E{0(t)} « 0 


E{Ut)C^(s)} - H6(t-s), E{0(t)0'^(s)} - 06(t-s) 


5 > 0 , 


0 > 0 


and 5(t), 0(t) are independent and Gaussian. The assumptions 

[MO, B(e)], (A(E), L(E)J controllable 

[A(e), C], [A(e), /q] observable 

are made, 0 ^ c ^ e^. 

As is well known, the -optimal control law is 
u(tj e) - - G(e)x(t> e) 

where 

-1 T 

G(e) - R B K(e) 
and K(e) satisfies 

0 - - K(e)A(e) - A'^(e)K(e) - Q(e) + K(e)B(e)R"^B'^{ e)k(e) 


The estLmate satisfies the equation 


(2.29) 


(2.30) 


(2.31) 


(2.32) 


(2.33) 


X(t» E) - A(E)x(t; E) ♦ H(E)(y(t> E) - Cx(t| E) ) + B(E)u(t» €) (2.34) 


where 

H(E) - E(E)C'0~^ (2. 

and £(e) satisfies 

0 » E(e)a^(e) + a(e)E(e) + l(e)5 l^(e) - E(e)c^r ^cE(e) (2. 

At this point, we ^u:e ready to apply Theorem 1 to approximate the 
controller (2.31), (2.34) by a two time scale controller. Because of 
Theorem 1, any system that has a stochastic degenerate system in ccxnmon 
with the optimal closed- loop system will be asympotically optimal. The 
optimal closed-loop system cem be vrritten 

r Jc^(t;E) ] 


ftj^(tj e) 

dc.(tf E) 
ex2(tiE) 








A^^-Bj^Gj^(e)-Hj^(e)c^ H^(e)C2 

A 21 -B 2 Gj^(e) A ^2 

eH_(e)c. A -B G (e)-eH.(e)c, eH_(e)c. 

21 ^'^21 21 22 


Aij-B^C^(C)-H^(E)Cj 

-BjGj(e) 

Aj2-BjGj(C)-eH2(E) 


Xj^(t; e) 

Xj^(t;e) 

X 2 (t;e) 

X 2 (t;e) 


0 H^(e) 


0 «H 2 (e) 


(2.37) 


Note that the stochastic degenerate system can be obtained by eliminating 
^ 2 ' ^2 equations for x^, x^ using the algebraic relations that 

result when C is set equal to zero in the left hand side of (2.37). Of 
course, the resulting system cannot be implemented since the value of e in 
the x^ and X 2 equations is not a design parameter. 

An implementable system that has a stochastic degenerate system in com 
mon with (2.37) is obtained as follows. Assume that (A^.,-B.^G^ (e)-H.,(e)C^)“^ 


22 2 2 


exists. Set e * 0 in the left heuid side of only the X 2 equations of (2.37) 
to obtain: 


X 2 - -(A 22 “ B2G2(£) - H2(e)C2) ( <A2^-B2G^ (e) - H2(E)Cj^)Xj^ + 

+ H2(E)Cj^Xj^ + 82(6)02X2 + 82(6)6] 

Substitute into the x^ equation to obtain 


^ljj(t; e) ■ ^iiD^ld^^’ 


(2.38) 


(2.39) 


where 


Aiid(6) - A^j^-B^G^(e) - H^(e)C^ - 

(A 22 -B 2 G 2 ( 6 )-H 2 ( 6 )C 2 )"^(A 2 i - B^G^(E)- 82 ( 6 ) 0 ^) (2.40) 


Hid(E) - H^(£)-(Aj^2"®1°2^^^“”i^^^^2^ ^'^22"®2®2^^^"”2^^^^2^' ” 2^^^ 


82 ( 6 ) - 6 H 2 ( 6 ) 


(2.41) 

(2.42) 


Baued on this analysis, the following suboptimal closed-loop system is 
obtained. 

ilD<t,E) 

Xu,<tlE) 

6 lt- (tie) 


6 ^ 


I 


11 


"21 

H2(E)C^ 


-BiG^(C) 

-B2G^(£) 

A2i-B2Gi(e)-H2(e)C^ 


12 


22 

H2<^)C2 




-B2G2(C) 






- "l 

0 


+ 

0 


Xjp(CE) 


^2 

0 



^ 0 

H^(C) -* 


-C(t)- 

-0(t). 


(2.43) 


Clearly, (2.43) is not the stochastic degenerate system of (2.37), but has 
been constructed to have a stochastic degenerate system in common with (2.37) 

Theorem 2 . 

Consider the LQG problem defined by (2.26) - (2.20). Assume that the 
controllai ility and observability assumptions (2.29), (2.30)hold and that 
(A^2 " B^G^ (e) -H2(eX^2^ exists, 0< e £ e^. Then the suboptimal closed-loop 

system (2.43) is asymptotically optimum in the following sense: 


lim 

e^O 


e{(x, (t»e) -X, (tjE) ) (x (t;£)-x (t;e) ) } » 0 

1 Id 1 id 


lim E{(x^(t;e)-Jc^j^(t|C)) ($c^(t;e)-Xj^jj(t;e))'^} - 0 

e^o 

lim E{(X2(t?e)-x (t;e)) (X2(t;e)-X2j^(t;e))^} » 0 

C*0 

lim E{(x2(t;e)-X2jj(t;e)) (X2(t;e)-X2jj(t;e))'^} = 0 
e^O 


Proof 

As noted above, (2.37) and (2.43) have a common stochastic degenerate 
system which governs their behavior as e 0. Therefore, all that is 
required is the verification of the hypotheses of Theorem 1. 

The required continuity properties follow from [18, 20]. Note that 
under the stated controllability and observ2d)ility assumptions the degener 
ate and boundary layer system controllability and observability conditions 
required are satisfied [24]. The required stability properties are es- 
tedslished as a consequence of the stability of the closed loop LQG design 
and a previously quoted result [22]. Full details can be found in [23]. 


Remar )^s 

1. Note that the above results are not the most general possible, since 
the time-varying and finite horizon case could probably also be solved. 
However, the aUx>ve results are of the greatest practical interst. 

2. Extension to more than two time scales in straightforward [22]. 
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3. Investigation of the rather ad hoc assumption (A 22 “B 2 (E)d 2 (E) (E) C^) 

exists for 0 < e ^ Eq would be of theoretical interest since verification 
is difficult. As a practical matter, this issue is less important, since 
invertibility for the value of e of interest is easily checked. The per- 
formance of the two time scale controller can then be directly assessed. 

4. The proposed suboptimal controller is illustrated in Figure 1. Notice 
that there is a unidirectional interface between the slow and fast filters. 
Thus there is opportunity for considerable reduction in on-line comput- 
ational effort since the two sets of filter equations can be numerically 
integrated in different time scales (i.e., with different step sizes). 

5. The above development assumed, for simpliciLy, that the optimal gain 
matrices H(e), G(e) were computed. In fact, by a more elaborate analysis, it 
is possible to show that only asymptotic approximations to H(E), G(E) near 

e “ 0 need to be computed [22]. It is possible to compute these approxi- 
mations by the methods extensively studied in [18-21] , thus obtaining a 
potential reduction in off-lint as well as on-line computation. 

6. Notice that it is the two time scale nature of the closed loop 
system that is required. In the above analysis, as £ -► 0 the closed loop 
system automatically has fast and slow nodes. In a physical system e has a 
fixed, non— zero value. Consequently, it is possible that the open loop 
dynamics can have fast and slow modes, but in the closed loop these slow 
modes au:e eliminated. 

TWO TIME SCALE CONTROL OF AIRCRAFT 
LONGITUDINAL DYNAMICS 

The particular problem to be addressed in this section is the design of 
a feedback control system for the longitudinal dyneunics of an F-8 aircraft. 
Specifically the controller must produce elevator commands to keep the air- 
craft in steady level flight in the face of wind disturbances. For simplici- 
ty, the wind disturbance is modelled as white. 

The equations of motion of an airplane are a set of coupled nonlinear 
equations in the longitudinal and lateral state variables. If the equations 
are linearized about nominal state and control variables, then the result- 
ing linear equations are found to approximately decouple into separate sets 
for the longitudinal and lateral dyneunics. See [1] for an excellent discus- 
sion of the modelling issues. 

The aircraft's longitudinal variables are 

r V 


L q 

where 

V » horizontal velocity deviation in ft/sec 

Y * flight path emgle in radians, 
a « angle of attack in radians, 

q ■ pitch rate in rad/sec. 

6 ■ is the elevator deflection in radians 

The Interpretation of these variables is given in Figure 2. Table 1 gives 
the system matrices. It is assumed that velocity and pitch rate measurements 
both corrupted by wideband noise, are available. 

Figures 3-4 show the system response to an initial pitch 6(0) ■ 1®, 

V(0) ■ 100 ft/sec in the eibsence of the wind disturbance. The two time scale 
behavior is well illustrated here. Table 2 gives the system's eigenvalues 
and eigenvectors. Note that the variables V, y dominate the slower phugoid 


, u - 


I 


mcxie, and the variables a,q dominate the faster short period mode. The 
physical nature of these oscillations is beautifully described on pages 
320-328 of (!]. 

From the above discussion and inspection of the system A matrix, the 
equation? for V and y are the logical candidates for the slow dynamics and 
the equations for a and q are suggested as the fast dynamics. For determi- 
nation of e, the following procedure is suggested. Note that the state 
equations can be written 



T T 

where *= tV yJ » ^2 ^ ^ ^ matrix 


GA 


22 


-.012 .0100 

-.0901 .0069 


has eigenvalues comparable to those of A,,. Determination of a value of e 
is actually not required for design of the two time scale controller, but 
is useful for judging the appropriateness of the approximation E = 0. 

Before proceeding with the regulator design, a few remarks are in 
order. First, the aircraft longitudinal variables more often include the 
pitch 9 instead of the flight path angle y. It was only after considerable 
difficulty that the above formulation, in which the choice of fast and slow 
variables is clear, was hit upon. Second, the extensive literature on 
singulau: perturbation theory contains almost no discussion of the choice 
of fast and slow variables, or the choice of e. But determination of these 
quantities is thr first problem one has to face up to in applications of 
the theory. 

A LQG controller was designed for the Q and R matrices Q ■ diag 
[.01,0,3260,3260], R • [3260J .using standard routines [29]. Design goals 
were to (i) achieve a deunping ratio ^ >. 707 for both modes (ii) reduce 
the state variables response to the wind disturbance. Closed-loop eigen- 
values and RMS state variable amd estimation error standard deviations are 
given in Table 3. As noted in remark 6 of Theorem 2, it is critical that 
the open loop separation of modes be present also in the closed loop. From 
Table 2, the optimal design has this property. The closed-loop eigenvalues 
of the two time scale controller euid corresponding RMS state variable and 
estimation error stemdard deviations are also listed in Tedsle 3. Note the 
generally good correspondance between the optimal and suboptimal designs. 


SUMMARY AND CONCLUSIONS 

This paper has considered the reduction in on-line computational ef- 
fort for an LQG design with fast and slow closed loop modes. Together with 
the results of [18-21], this paper demonstrates that the singular pertur- 
bation approach to the LQG problem offers the potential of near optimal 
performance with reduction in both on-and off-line computation. 

The design procedure of this paper has been applied to control of the 
longitudinal dynamics of a jet aircraft. A two time scale design was 
obtained with performance extremely close to that of the optimal design. 
Note that recent proposals for adaptive flight control systems require 
multiple Kalman filters running in parallel [30, 31], so that reduction of 
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on-line computation is of definite interest. 

Several directions for future research are evident. First, procedures 
for systematically picking the fast and slow variables of a system with fast 
and slow modes, as well as for determining explicitly f, would be highly 
useful in applications. Second, note that the design of Figure 1 has an 
interesting hierarchical structure. In fact, as pointed out in [14], if a 
system is composed of a number of fast subsystems with a slow interccnnect- 
ing equation, a decentralized, hierarchical design is naturally obtained. 
Therefore, the results of this paper are of potential interest in hierarchi- 
cal systems theory. Finally, note that singular perturbation theory is only 
one approach to the multiple time scale phenomenon. The results of Ramnath 
[32], for example, provide another approach that could be exploited in 
control theory. 
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APPENDIX 

Proof of Theorem 1: 

From (2.7), (2.14), (2.15) the following equations are obtained 


x^(t»G) 

ex^ (t;£) 

Xj^(t;£) 

€^2 



0 


A^j(£) 


*22d<^> 

0 

0 

0 

0 


Ajj(E) 

0 

0 

Ajl(e) 

AjjCO 


Xj^(t>e) 
X2(t;C) 
Xj^ (t;G) 
X2(t;G) 


where 


Lj^(e) 


L^(G) 

L2(g) 


C(t) 


(A.l) 


Ai3(e) - A^2<^J^22^^^^^21^^^ 


Li(e) - (£)L 2 (e) 


(A. 2) 
(A. 3) 


Define 


u(tie) 


£33(t;C) 




Ej3(t,e) 

Z2^(t;e) 

3^<t,e) 

E3j(t,e) 

l33<t,e) 

Z3^(t,E) 


1 ^ 2 (t,c) 

Z^3(t.e) 

5:44 


and let 


x^ (t;E) 
X 2 (t;e) 
Xj^(t;e) 

X2(t;e) 


x^(t;e)X2(t;e)x^(t;e)X2(t,*e) 


(A. 4) 


'44 


(t;£) - £Z^^(t;£) 


(A. 5) 


12 


I 


Then the variance equations corresponding to (A.l) can be written 




+ Lj^(e)HL[(e) 


(A. 6) 




L^(e)HL^(e) 


(A. 7) 


^33 “ '*■ ^12^^^^43 ^ ^33^11^^^ ^34^12^^^ L^(e)cLj^(e) 

(A. 8) 

^^12 ” ^lld^^^^l2 ^13^^^^32 * ^^12^^^^42 ^ ^11^21d^^^ ^12^22d^^^ 

(A. 9) 

^^14 “ '*' ^^13^^^^34 ■*■ ^12^^^^44 ^13^21^^^ ■*■ ^14^22^^^ 


+ L^(£)5L2(e) 


(A. 10) 


^^22 “ ^ 21 d ^^^^ l 2 ^ 22 d ^^^^22 ^ 21 ^ 21 d ^^^ ^ 22 ^ 22 d ^^^ 


(A. 11) 


^^23 ' ^21d^^^^l3 ^22d^^^^23 ^ ^^23^11 ^^24^12 


(A. 12) 


^^24 “ ^21d^^^^l4 ^22d^^^^24 ^23^24^^^ ^24^22^^^ 


(A. 13) 


€^34 »= GAj^^Ce)!^^ + ^12^^^^44 * ^33^21^^^ ^34^22^^^ '*’ Lj_(E)=L2^^^ 


(A. 14) 


^^44 • ^ 21^^^^34 ^ '' 22<^^^44 ^ ^^ 43 "' 21 ^^^ * ^44"^^^^^ L 2 (£)HL ^(£) 

(A. 15) 

To apply the theorem of Hoppensteadt to (A. 6) - (A. 15), conditions 

(The conditions are conveniently listed in [18]) must be verified. 

Conditions H.., H. and H are satisfied since (A. 6) - (A. 15) are linear 
11 IV V 

time-invaricuit differential equations with coefficients continuous in £. 

Condition is verified by noting that (A. 9) - (A. 15) have a unique 

solution, for e « 0 when ^13' ^33 given, as follows. Since A22(0), 

^22ci^^ * 1 ^® stable, ^23 '-’niquely determined by 


I 


I 


0 = + Lj(0)=Lj(0) 

0 - ^12*22d<°> 


(A. 16) 


(A. 17) 


° “^21(^°^^13^ ^22d^°^^23 


(A. 18) 


Therefore, Z,,, and I,, are the unique solutions of 

14 22 34 


0 = A^ 2 ( 0 )i^^ + 

° '''21d<°>!:i2*''22d«''>S2 * =^2l'‘nd*‘» * ^22''L"» 

0 - A^j( 0 )E^^ * * ^ 34 '' 22 “'' * 


(A. 19) 


(A. 20) 


(A. 21) 


Finally, 2^24 uniquely determined by 

° " *21d«»=^l4^''22d<'»=:24*^23''L'°>*^2/l2<‘» 


(A. 22) 


To verify H . , it must be shown that the solution defined by (A. 16) - 
(A. 22) is an asymptotically stable equilibrium of the following boundary 
layer system associated with (A. 6 ) - (A. 15). 


^ ^12 = ^12''Ld<°> 


(A. 23) 


A = a^2(0)Z,, * Zi3A^,(0) * E3,A^2(0) 4 E^(C)=L^(0) 


(A. 24) 


H ^22 “^ 21 d ^°^^12 ■*‘^ 22 d ^°^^22 ^ 21 ^ 2 ^^°^"^ ^ 22 ^ 2 d 


(A. 25) 


dT ^23 “^21d^°^^13 ■^^22d^°^^23 


(A. 26) 


^ ^24 “^ld^°^^14 ^W°>^24 ^ ^23^L<°> * ^24"L^°> 


(A. 27) 


^^34 = ^2<°)^44 " " I-i(0 )EL^(0) 


(A. 28) 


^ ^44 * '‘22«»^44 * ^44^22 


(A. 29) 


The stability of (A. 2 3) - (A. 29) is readily verified using the stability 
of A 22 ( 0 )» ^ 22 d^^^ wor)cing through the equations in the order indicated 

in (A. 16)- (A. 22) . 

The final conditions deal with the degenerate system of (A. 6 ) - (A. 15). 
Putting e ■ 0 in (A. 6 ) - (A. I 5 ) , 6 md after som« algebra using (A. 16) - (A. 19) , 

(A. 21), the following equations are obtained: 
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E 33 - 

+ L^(0)E^L^(0) 















(A. 

30 ) 

'"il 

( 0 )- 

■* 12 <*» 

*^ 2 ^ 

(0)A2i 

( 0 ) 

’*• 31*^31 

**11 

( 0 )- 

0 

CM 

< 

t 

*22 

( 0 ) 

^21 

( 0 ) 

)" 














(A. 

31 ) 

**11 

( 0 )' 

-* 12 *°> 

* 2 ^ 

<°'*21 

( 0 ) 

'=^ 11*^11 

(All 

( 0 )- 

- A ^ 2 < 0 ) 

*- 2 ^ 

( 0 ) 

^21 

( 0 ) 

)^ 














(A. 

32 ) 


Since A (0) - A (0)A ^{0)A (0) is stable, H. and H . are satisfied. 

XJL Xm 2«2X 1 V i 

Thus by Hoppensteadt ' s Theorem, lim Z j^(t;C) and lim l 22 (t;e) are 

e-K) e-^0 

equal to the solutions of (A. 32) and (A. 20). Since 

x^(0,E) = Xj^(O) - x^j^(O) * 0 (A. 33) 

by choice of the initial conditions of the stochastic degenerate system, it 
follows that the initial conditions of (A. 32) are 


Z^^(0; 0) = 0 


(A. 34) 


Thus 


lim Z (t»e) = 0 (A. 35) 

From (A. 35), (A. 17) and (A.20) , 

lim Z 22 (t;E) = 0 (A. 36) 

C-+-0 

Note that (A. 30) is valid at t « 0 by choice of initial conditions. 
Equations (A. 34) and (A. 36) are the desire 1 result. 
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eigenvector conponents 



f 


A = 


B = 


C * 


-1. 357x10 
1.200X10 
-1.212X10 
5.700x10 


-3.220X10^ 

0.000 

0.000 

0.000 


-4.630X10^ 

1.214 

-1.214 

-9.010 


0.000 

0.000 

1.000 

-6.696x10 


- 4 . 330 x 10 "^ 


-4.630X10^ 

1 . 394 x 10 ”^ 


1.214 

- 1 . 394 x 10 "^ 

L = 

-1.214 

-1. 577^10"^ 


-9.010 





5 “ [3.150X10*'^] 


_ 



- 


-4 


0.0 

0 

• 

0 

0.0 

1.0 

0 = 

6.859X10 

0.000 

1 

1.0 

0.0 

0 

0 

0.0 

id 

0.000 

4.000X10 


Table 1. System Matrices 


eigenvalue 


-.94 ± j.2.98 -.0075 ± j.076 


1.3 ± j.4.7 -4.3X10^ 4 jl.lxiO^ 


• 5 x 10 "'^ + j.l3 


-.15 ± jl.O 


a 


.33 ± j.lO 


-.02 T j.005 


-.21 ± jl.O 


-.00 j.02 


Tcjble 2. Open Loop Eigenvalues 
and Eigenvectors. 







-.94 ± j2.98 -.075 ± jO.76 
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